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In this communication we investigate the quantum statistics of three harmonic 
oscillators mutually interacting with each other considering the modes are initially in 
Fock states. After solving the equations of motion, the squeezing phenomenon, sub- 
Poissonian statistics and quasiprobability functions are discussed. We demonstrate 
that the interaction is able to produce squeezing of different types. We show also 
that certain types of Fock states can evolve in this interaction into thermal state and 
squeezed thermal state governed by the interaction parameters. 
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I. INTRODUCTION 

In the field of nonlinear optics there are two basic parametric processes, namely, para- 
metric amplifier and frequency converter. The parametric amplifier is designed to amplify 
an oscillating signal by means of a particular coupling of the mode to the second mode of 
oscillation, the idler mode. The coupling parameter is made to oscillate with time in a way 
which gives rise to a steady increase of the energy in both the signal and idler modes 
The parametric amplifier is the source of squeezed light (i.e. the states of light with reduced 
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fluctuations in one quadrature below the level associated with the vacuum state [2^). To 
be more specific the degenerate and non-degenerate parametric amplifier are the sources for 
single-mode and two-mode squeezing, respectively. The lossless Hamiltonian describing such 



a system (non-degenerate case) is given by 
H ^ 

^^Wja^fij + i\i {aia2 exp[i(a;i + uj2)t\ — h.c.} 



n 



where aj and hj are the creation and annihilation operators, Ai is the coupling constant 
including the amplitude of the pump, and Wj, j = 1, 2 are the field frequencies. On the other 
hand, the parametric frequency converter is described by a process of exchanging photons 
between two optical fields of different frequencies. The Hamiltonian representing such a 



system reads 



h 



jdj + 1X2 \ a\a2 exp[i(cJi - ^2)^] - h.c. \ . 



;i-2) 



This model can be applied to describe various optical phenomena, e.g. to find analogies 



between frequency converter and beam splitter 



electromagnetic field {4], and Raman scattering [3|,|5|. The quantum properties of parametric 
frequency converter are discussed in 6|. Furthermore, some authors studied this model as 
the lossless linear coupler, e.g. in j?]. In this situation the model is considered to be 
represented by two electromagnetic waves which are guided inside the structure consisting 
of two adjacent and parallel waveguides; the linear exchange of energy between these two 
waveguides is established via the evanescent field [8|. 

The derivation of the above Hamiltonians can be obtained from the quantization of the 
cavity modes in a volume V . In this case the total energy of the field is given by js, 9| 

H =h Uty(e£-+'^')d'V, (l-3a) 
where £ and 1-L are the electric and magnetic fields, respectively, and e is the dielectric 
constant. To provide coupling between the various cavity modes we have to make the 
dielectric constant e varied according to 

e(i?, t) = 1 + A tf{R) cos{LOit + 0,), (1.36) 
where Ae <^ 1 and f{R) is a function of the position vector i?, while 0j is an arbitrary pump 
phase. From the relations (1.3a) and (1.3b), and in terms of the boson operators, one can 
write the Hamiltonian (1.3a) in the following form 



r 



two-level atom driven by a single mode of 



H 



Y^iUJi{a^a + I) - Y^iim^im cos{ujit + 0i)(a[ - a/) (a. 



where kim is the coupling coefficient having the form 



7t _ 



;i.3c) 
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km = j^iuiUm)^ f^^^.^y f{r)ui{r)um{r)dV, {1.3d) 
where the function ui{r) represents the normal mode and satisfies 

VA VAm^W = (^)V(r). ^ (1.3e) 

It is convenient to refer to the recent work of Garrett et al. |10i], who reported an 
experimental demonstration of phonon squeezing in a macroscopic system by exciting a 
crystal KTqOs with an ultrafast pulse of light. The Hamiltonian relevant to their work can 
be adjusted to be essentially equation (1.3c). Thus if one chooses the function f{R) so that 
kim 7^ 0, this will leave an infinite number of modes coupled. Hence by a proper choice of 
the pump frequency Ui, the interacting modes will be limited to two modes, and then under 
certain conditions one can obtain the Hamiltonians (1.1) or (1.2). 

The optical processes involving the competition between frequency converter and para- 
metric amplifier are of interest from theoretical and experimental points of view, which can 

12 1 and nonlinear coupler {13!. 11^. where the nonlinear 



be seen in three-mode interaction 
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interaction couples different photon modes and leads to the energy transfer between modes. 

Photons in some modes may be annihilated, while those in the other modes created, and 

hence the photon distribution is disturbed. In each time the rate of energy transfer between 

the modes depends on the statistical properties of the light fields. This in fact encouraged 

us to follow up the procedure given above and adjust the pump frequency to extend the 

number of modes to be three instead of two. This can lead to the following Hamiltonian 
H ^ 

— = [^Wja]dj] + i\i {0,102 exp[i(wi + ^2)^] - h.c.} 

+1X2 jaicig exp[i(a;i — uj3)t] — h.c.j + iX^ 1 0,20-3 exp[i(a;2 — W3)t] — h.c.j , (1.4) 
which represents a model of three boson field modes, designated by the annihilation opera- 
tors dj,j = 1,2,3, and an intense monochromatic light wave inducing parametric coupling 
between them. In this Hamiltonian Xj are the coupling constants including the pump am- 
plitude and proportional to the second-order susceptibility of the medium x^"^^ are the 
natural frequencies of oscillation of the uncoupled modes and h.c. is the Hermitian conju- 
gate. This interaction can also be established, e.g. by means of a bulk nonlinear crystal 
exhibiting the second-order nonlinear properties in which three dynamical modes of frequen- 
cies Ui,u)2,u)3 are induced by three beams from lasers of these frequencies. When pumping 
this crystal by means of the corresponding strong coherent pump beams, as indicated in 
the Hamiltonian, we can approximately fulfil the phase- matching conditions for the corre- 
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spending processes, in particular if the frequencies are close each other (biaxial crystals may 
be helpful in such an arrangement). Also a possible use of quasi-phase matching may help 
in the realization, which is, however, more difficult technologically 15|. There is another 
possibility to realize the interaction described by the Hamiltonian (1.4) using a nonlinear 
asymmetric directional coupler as shown schematically in Fig. 1. 

For completeness, Hamiltonian (1.4) is a generalization for several models used in the 



literature, e.g. [l|, [ll 



16| . and it is invariant under the transformation cli < — > a2 with 



A2 i — A3. So that the equations of the 2nd mode can be obtained from those of the first 
one. 

a2(0) 



am 



am 



X^'^ Ai 
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.ai( 



X 



(1) 



Fig.l Scheme of realization of interaction in (1-4) using a nonlinear asymmetric direc- 
tional coupler which is composed of two optical waveguides fabricated from first-order (x^^'^ ) 
and second-order (x^'^^ ) materials, where x designate susceptibility. Modes 1 and 2 propagate 
in the first waveguide and mode 3 in the second waveguide. The interaction between modes 1 
and 2 is established by strong pump coherent light, which is not indicated in figure, with the 
coupling constant Ai. The interactions between mode 3 and modes 1 and 2 are established 
linearly with the coupling constants A2 and A3, respectively. The beams are described by the 
photon annihilation operators as indicated; z = vt is the interaction length and we assume 
that all beams have the same velocity v and the length of the waveguides is L. Outgoing 
fields are examined as single or compound modes by means of homodyne detection to observe 
squeezing of vacuum fluctuations, or by means of a set of photodetectors to measure photon 
correlations, photon antibunching and sub-Poissonian photon statistics in the standard ways. 
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In this article we investigate the statistical properties for the different modes controlled 
by the Hamiltonian (1.4). We consider the modes are initially prepared in Fock (number) 
states where these states represent the most basic quantum states and are maximally distant 
from what one would call a classical field. More illustratively, they can always exhibit sub- 
Poissonian statistics, however, they provide no information on the phase since the number 
of photons is quite certain. Also any states can be expressed in terms of these states with 
corresponding weighting distributions. Further, there are various proposals how to generate 
such states. For instance, these states can be prepared in a nondemolition measurements 



3 



17-11 or m the micromaser 



20| . During this work we use a finite number of quanta in each 



modes, which is relevant to the experimental realization. Also we treat the model as three 
harmonic oscillators interacting in a nonlinear crystal in a sense of a competition between 
two frequency converters and one parametric amplifier, as indicated in the Hamiltonian 
(1.4). Under these circumstances the model is motivated by various interesting results. For 
example, we show that this model can be used as an effective source for squeezed light addi- 
tionally to those systems given in the literature earlier 21|. Also we show that it can yield 
thermal light and squeezed thermal light. Moreover, we demonstrate how the qualitative 
behaviour of the single nonlinear process such as parametric amplifier or frequency converter 
can be substantially changed in the presence of the other nonlinear processes. 

The plan of the work is as follows: In section 2 we introduce the solution of the equations 
of motion, section 3 is devoted to a discussion of nonclassical phenomena such as quadrature 
squeezing, second-order correlation function and the violation of Cauchy-Schwarz inequality, 
in section 4 we demonstrate the quasidistribution functions, and finally we summarize main 
conclusions in section 5. 



II. SOLUTION OF THE EQUATIONS OF MOTION 

In this section we shall introduce the exact solution for the equations of motion in the 

Heisenberg picture for the Hamiltonian (1.4). 

The dynamics of the system is described by the Heisenberg equations of motion which 

for any operator 6 are given by 
dO dO 1 

where [...,...] represents the commutator. 
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Therefore the propagation of the field operators, which are defined in the slowly varying 
forms {dj — Aj ex.p{—iujjt), j — 1, 2, 3), can be written in the following form: 



~df 

dA2 

dt 
dA^ 

dt 



= -\iA\ - Aaia, 
= -Alii - Aaia, 



A2A1 + As^. 



(2.2a) 
(2.26) 
(2.2c) 



These are three equations with their Hermitian conjugates forming a closed system which 
can be solved easily as 

Mt) = «i(o)/i(i) + d\{Q)h{t) - dMhit) 

-dl{o)U{t) - d,{o)Mt) - dl{o)Mt), (2.3a) 



Mt) = d2{0)9i{t) + 4(0)^72(t) - d,{0)gs{t) 
-dl{0)g4{t) - ds{0)g,{t) - a^(0)^6(i), 



(2.36) 



Mt) = o,3{0)hi{t) + dl{0)h2{t) + d2{0)hs{t) 

+dl{0)hi{t) + ai(0)/i5(t) + d\{0)he{t), (2.3c) 

where the exact forms for the time-dependent coefficients, i.e. for fj{t),gj{t),hj(t), which 

include all information about the structure of the model are complicated for the general 

case, we give their form only for the special case when A2 = A3 = A (hence fj{t) = gj{t)) , 

which will be frequently used here: 
1 



/2,4(i) 

hit) = 



^2 



A 



cosh(^^) cos{kt) + ^ sinh(^^) sm{kt) ± cosh(Aii) 

sinh(^^) cos(kt) + ^ cosh(^i^) sm( kt) ^ sinh(Aii) 
2 2A; 2 



.Alt, 



^ cosh(^)sin(A;t), fe{t) 



A . Alt . - 
Y sinh(-— ) sm[kt), 



(2.4) 



hi(t) = cosh(^^) cos{kt) — ^ sinh(-;^) sin(A;t), 

A t X t 

h2{t) — — sinh(^-) cos(^t) + ^ cosh(-^) sin(^t), 

h3{t) = /i5(t) = ^ cosh(^^) sin(A;t), /i4(t) = he{t) = — ^ sinh(^^) sin(A;t), 

4 • 



(2.5) 



where k — \/2X^ — ^. In fact, the nature of the solution can show how the interaction 
does work. That is the time-dependent coefficients contain both trigonometric and hyper- 
bolic functions. Consequently, the energy associated with the propagating beams inside the 
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nonlinear crystal can be amplified as well as switched between modes in the course of time. 

On the basis of the well known commutation rules for boson operators, the following 
relations can be proved between the time-dependent coefficients 

flit) - fm + m - m + m - m = i, 

mg.it) - h{t)g,{t) + h{t)g2{t) - U{t)g,it) - Ut)g,{t) + mg,{t) = 0, 
mgsit) - f2{t)g,{t) + h{t)g,{t) - Ut)g,it) - h{t)g,{t) + h{t)g,{t) = 0. (2.6) 
The remain relations can be obtained from (2.6) by means of the following transformations 

(/l(t),/2(t),/3(t),/4(t),/5(t),/6(t)) 

< — > [-93{t), -giit), -gi{t), -g2it),g5{t),gG{t)) 

M (^h{t),he{t),-h{t),-K{t),-h{t),-h2{t)y (2.7) 
Based on the results of the present section, we can study the quantum properties of the 
evolution of the different modes in the model when they are initially prepared in Fock states. 



III. NONCLASSICAL PHENOMENA 



Since the photons produced in a nonlinear optical process are known to possess unusual 
correlation properties resulting in many of nonclassical aspects of the radiation field. There- 
fore, we investigate some of these properties for the model under consideration in the present 
section. Our attention is focused on squeezing phenomenon, sub-Poissonian photon statistics 
and the violation of Cauchy-Schwarz inequality. 



A. Squeezing phenomenon 



The concept of squeezing of a quantum electromagnetic field has been given a great deal 



ne 

3 



of interest in the view of the possibility of reducing the noise of an optical signal below t 
vacuum limit and the possible potential application in optical communication networks 
and gravitational wave detection 22|. This light can be measured by homodyne detection 
where the signal is superimposed on a strong coherent beam of the local oscillator. The 



23|. 



generation of such light has been reported in various schemes 

Generally there are several definitions for squeezing, e.g. see [2l|]. Here we investigate 
some of these types for the model under discussion, in particular, single-mode, two-mode and 
three-mode squeezing as well as sum-squeezing 2J] when the modes are prepared initially in 
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Fock states. The starting point for this study are the two quadratures X and Y which are 
related to the conjugate electric and magnetic field operators E and H. They are defined 
in the standard way. Assuming that these two quadrature operators satisfy the following 
commutation relation 



X,Y 



C, (3.1a) 



where C may be an operator or a c-number with respect to which kind of squeezing we 
consider, the following uncertainty relation holds 

((AX)2)((AFf)>M^, (3.15) 

where ((AX)^) = (X^) — (X)^ is the variance. Therefore, we can say that the model 
possesses X-quadrature squeezing if the S'-factor [251] . 

Sit) = ((^Mm^oMm 

0.5|(C)| 

satisfies the inequality — 1 < 5 < 0. Similar expression for the F-quadrature (Q-factor) can 
be obtained. 

Firstly, for three-mode squeezing we define the two quadratures as 

Mt) = llAiit) + A^it) + As{t) + A\{t) + Al{t) + Al{t)], (3.2a) 

%{t) = l[ii(t) + A,it) + A,{t) - Alit) - Alit) - Alit)], (3.26) 

where the subscript 3 in the left-hand side stands for three mode case. Using (2.3) together 
with the definition of the variance, for initial Fock states we arrive at 

4((AX3(t))2) = (1 + 2ni)[/i(t) + /2(t) + hit) + Kit) - g,it) - g,it)]' 

+(1 + 2n2)[/3(t) + /4(t) - giit) - g2it) - hit) - hit)f 

+ (1 + 2n3)[/5(t) + m + g,it) + g^it) - h{t) - h{t)]\ (3.3a) 



mA%it)f) = (1 + 2ni)[/i(t) + g,it) + hit) - hit) - g,it) - h{t)f 
+ (1 + 2n2)[hit) + g2it) + hit) - f,it) - (71 (t) - h{t)f 

+ (1 + 2n,)[hit) + g,it) + h{t) - hit) - g^it) - h{t)f-. (3.36) 
fij is the mean photon number for the jth mode. The expressions for the single-mode 
and two-mode squeezing can be obtained easily from (3.3) for two quadratures defined in 
a manner analogous to (3.2) by dropping the coefficients of the absent mode, e.g. for the 
1st mode, single-mode squeezing can be obtained by setting gjit) = hjit) = 0,j = 1,2,3 in 
(3.3). It should be taken into account that (7 is a c-number in this case and equals |, 1, | 
corresponding to the single-mode, two-mode and three-mode squeezing, respectively. 
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We start our discussion when A2 = A3 = A. In this case the quadrature variances for the 
3rd mode (single-mode squeezing) are 

((AXi(t))2) = 1 |2(1 + n,+ n2)[^ sm{kt)f 

+ (1 + 2n3)[cos(H) + ^sin(H)]H exp(-Ait), (3.4a) 

2k J 

((Afi(t))2) = i|2(l + ni + n2)[^sin(H)]2 

+ (l + 2n3)[cos(A;t) - ^ sin(H)]H exp(Ait). (3.46) 

2k J 

The immediate conclusion can be drawn from these expressions which is that there is squeez- 
ing in the X-quadrature after switching on of the interaction by a suitable time as a result 
of the fact that the Fock states are not minimum-uncertainty states. Indeed, this is an 
interesting result because most of the nonlinear optical processes including correlation of 
modes are unable to provide such a property. Furthermore, under this condition similar 
expressions, i.e. (...) exp(— Ait) for X-quadrature and (...) exp(Ait) for F-quadrature, can 
be obtained for the two-mode squeezing, in particular, between the 1st and 2nd modes and 
also for the three-mode squeezing. Also, this situation is still valid if one considers that 
the modes are initially in the states which are represented by a density matrix diagonal in 
the number state basis (thermofield states); in this case fij represents the average photon 
number for the jth thermal mode and therefore the behaviour of the model under considera- 
tion is close to that of a microwave Josephson-junction parametric amplifier [26|, in which a 
thermal input field has been introduced to the squeezing device and the generated field has 
exhibited noise reduction. Indeed, thermal noise is inevitable and hard to quench, so that it 
is more realistic to consider a thermal state instead of a vacuum state input to a squeezing 



device 



27|. More detai 



states can be found in 



s about the properties of squeezed number and squeezed thermal 
. We proceed by plotting the squeezing factors associated 



27 



with the X-quadrature, i.e. S'-factor, for the three types when the modes are initially in 
the Fock states as shown in Fig. 2, for the given values of the parameters. In this figure we 
have considered the single-mode (star-centered curve-3rd mode), two-mode (short-dashed 
curve-(l,2) modes) and three-mode squeezing (long-dashed curve). For the sake of compar- 
ison we adopted solid curve for squeezing factor of well-known squeezed number states 29| , 
i.e. for Si{t) = (1 + 2n) exp(— Ait) — 1, r = Ait is the usual squeeze parameter. From this 
figure we can see that for all cases at the first moments of the evolution there is no squeezing 
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owing to the nature of the Fock states. Further, the three types have approximately similar 
behaviours in a sense that they oscillate around the value of the usual squeezed number 
states, i.e. around the solid curve, as a result of switching of the energy between the inter- 
acting modes in the system. The significant remark is that they arrive to the steady state 
(maximum squeezing) at the large values of interaction times. 

From the above analysis it is clear that the 3rd mode possesses more pronounced nonclas- 
sical behaviour compared with the 1st and 2nd modes. This is related to the structure of the 
Hamiltonian (1.4) and can be interpreted as follows: modes 1 and 2 undergo amplification 
from the first process (i.e. from the parametric amplifier) and then the energy transfers from 
these two modes to the 3rd mode via the converter processes. Indeed, such a structure is 
the source of the single-mode squeezing which has not been seen before for the most of the 
important models involving three-mode interaction process, e.g. jlll. l3o|. 




time 

FIG. 2: Squeezing factor, 5„(t), for initial number states with mean photon numbers fij = 1, j = 
1, 2, 3, against the time for the single-mode (star-centered curve-3rd mode), two-mode (short-dashed 
curve-(l,2) modes) and three-mode squeezing (long-dashed curve) for (Ai, A2, A3) = (0.25,0.3,0.3); 
the solid curve is for squeezed number states, i.e. for Si(t) = (l-|-2n) exp(— Ait) — 1, Ait represents 
the usual squeeze parameter r. 



Secondly, we discuss sum-squeezing [24] for the model under discussion. Sum-squeezing 
effect is both higher-order and multimode phenomenon and the states join to this class are 
nonclassical states since its Glauber-Sudarshan P-function always does not exist and it is not 
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well defined function. Also, there is a connection between sum-squeezing and sum-frequency 
generation and consequently it can be converted into normal single-mode squeezing by an 
appropriate nonlinear optical process and this may be used for the detection. The importance 
of this type is related to the fact that squeezing can exist for the correlated modes even if 
the individual modes are not squeezed in the normal sense. 

Now the quadratures corresponding to the real and imaginary parts for sum-squeezing 



are given by 



24| 



n(t) = i[i,(t)i,(t)-i|(t)il(t)], 



(3.5a) 
(3.56) 

where j ^ k and the subscript s stands for sum-squeezing. In this case C = A^j{t)Aj{t) + 
Al{t)Ak{t) + 1, i.e. it is a state dependent. One should notice that the amplitude squared 



squeezmg 



31 1 is the degenerate limit of sum-squeezing, and also the operators Xs{t), Ys{t), C 



form a representation of the su{l, 1) Lie algebra and thus they have been used in the inter- 



ferometric measurements 



32|. 




-0.7 



time 



FIG. 3: Sum-squeezing factor for initial number states with mean photon numbers fij = 1, j = 
1,2,3, against the time for Xj = 0.1, j = 1,2,3 and for the single parametric amplifier (Y- 
component, solid curve), between modes 1,2 (Y-component, short-dashed curve) and between 
modes 1,3 (X-component, long-dashed curve). 



The expressions for S{t) and Q{t) for sum-squeezing are rather lengthy in this case, 
further, they can be calculated straightforwardly and are not very illuminating, so we shall 
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not quote their explicit forms. As is known the single parametric amplifier is the source of 
sum-squeezed light jssi, however, single frequency converter cannot produce such a light. 
The system under consideration can produce sum-squeezed light which is more effective than 
that of the single parametric amplifier and this refiects the role of competition of different 
processes in the model. These facts are clear in Fig. 3 where the sum-squeezing factor is 
plotted against the time for the shown values of the parameters. Solid, short-dashed and 
long-dashed curves show that sum-squeezing factor corresponds to the F-component of single 
parametric amplifier, the y-component of (1,2) compound mode and the X-component of 
(1,3) compound mode, respectively. From this figure one can observe that for the single 
parametric amplifier (solid curve), it starts from squeezing bound, goes to its maximum 
squeezing value (~ 12%) and then squeezing values decrease and eventually comes back to 
the initial values for large interaction times. However, for (1,2) mode, squeezing appears 
after switching on of the interaction as before, then its value increases in the course of time 
and eventually locks to give 50% squeezing in the large time of interaction domain. It is 
important to point out neither the 1st mode nor the 2nd one is squeezed in the normal sense 
in this case. We proceed by discussing the quantity for (1,3) mode (long-dashed curve). 
One can observe that squeezing is periodically established displaying its maximum value 
in the large interaction time and then squeezing disappears (this is not indicated in the 
figure). The maximum value of the sum-squeezing in this case is more pronounced than in 
the previous ones, in particular, up to 70% squeezing can be obtained. We should remind 
here that the 3rd mode can produce squeezing in the normal sense (cf. (3.4)). Consequently 
we can conclude that there is no connection between sum-squeezing and normal squeezing 
for the models including correlation between modes 2J]. It is important to point out that 
the behaviour of the sum-squeezing in (2, 3) mode is similar to that for (1, 3) mode and 
this can be recognized from the Hamiltonian (1.4). Final remark, the squeezing values are 
sensitive to the initial mean photon-numbers of the modes which can vanish completely for 
its large values. 



B. Sub-Poissonian statistics 



Sub-Poissonian light is an example of nonclassical light which can be measured by pho- 
todetectors. A state (of a single mode for convenience) which displays sub-Poisson statistics 
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is characterized by the fact that the variance of the photon number ((Anj(t))^) is less than 
the average photon number {hj{t)) = {Aj(t)Aj(t)) . This can be expressed by means of the 
normahzed normal second-order correlation function as 



1 + 



(3.6) 



('2') 

where the subscript j relates to the jth mode. Then it holds that g- (t) < 1 for sub- 



Poissonian distribution, gf\t) > 1 for super- Poissonian distribution and when g)l^\t) = 1 
Poisson distribution of coherent photons occurs. Also the system exhibits maximum sub- 

(2) 

Poissonian statistics when g^ (t) = 0. Here we may mention that Fock state, chaotic field 
and coherent state are representing good examples for sub- Poissonian, super-Poissonian and 
Poissonian statistics, respectively. Moreover, the generation of sub- Poissonian light has been 
established in a semiconductor laser 



in the microscopic regime 



34j and in the microwave region using masers operating 



35| (for review see 



On the other hand, it is worth referring to photon antibunching (bunching) phenomenon 
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which is generally not equivalent to sub- Poissonian (super-Poissonian) photon statistics 
38| . The basic formula to_ study this phenomenon is the two-time normalized intensity 
correlation function 



38 



39|, which can be represented as the joint detection probability 



of two photons, one at time t and another at time t + t, where r is the time interval 
between the two-photon detection process. Nevertheless, the bunching/antibunching and 
super-/sub-Poissonian statistics are identical to the stationary fields. For more details about 
the connection between these two phenomena, different definitions as well as applications 
the reader can consult |39|. 



Now we investigate the sub-Poissonian statistics for our model when all modes are initially 
in the Fock states as we did for squeezing. In this case the moments {A[{t)Ai{t)) and 
{A\^{t)Aj{t)) for the 1st mode read 

{A\it)Mt)) = [fm + fm]ni + [fm + fm]n2 

+[mt) + fMns + flit) + m + (3.7a) 



{Af{t)Al{t)) = ft{t)n,{n, - 1) + f^{t){n, + l){n, + 2) + ft{t)n,{n2 - 1) 
+ftit){n2 + l)(n2 + 2) + /5'(t)n3(n3 - 1) + ftmn^ + !)(% + 2) 
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+f,{t)f,mn, + If + [/i(t)/2(t)(2ni + 1) + /3(t)/4(t)(2n2 + 1)]^ 

+^mt)ni[f!{t)n2 + f!mn2 + l)] 

+4/|(t)[/3 W^2(rii + 1) + f!mni + 1)(^2 + 1)] 

+2(2^3 + l)/5(t)/6(t)[/i(t)/2(t)(2ni + 1) + Mt)Uit){2n2 + 1)] 

+4[/2(t)n3 + fim^3 + + mini + l) + + f.mn, + l)]; (3.76) 

as we have mentioned earlier nj is the mean photon number of the jth mode. The expressions 
associated with the 2nd and 3rd modes can be obtained from (3.7) using the transformations 
(2.7). 

Let us start our discussion when ni = n2 = = n and also consider the case of negligible 
contribution from parametric amplifier (Ai ~ 0). In this case the normalized normal second- 
order correlation function for the 1st mode reduces to 

g?\t) = 2 - (1 + i)[/f (t) + m + fM- (3.8) 

From the identity (2.6) one can show that < {ff{t) + /|(t) + /|(t)) < 1 and thus 

2 >g^;'\t) >{!--), (3.9) 
n 

which means that the model, in this case, can exhibit only sub-Poissonian statistics or partial 
coherence. 

Now we can investigate the evolution of Qj (t) when the modes are initially prepared in 
the number state |1, 1, 1) for the 1st and 3rd modes, as shown in Figs. 4a,b, respectively, 
for shown values. In Fig. 4a one can observe that g^\t) has smoothed behaviour (see, 
solid and long-dashed curves), i.e. it goes rapidly from values corresponding to the sub- 
Poissonian statistics to those for super-Poissonian ones and becomes almost stable for large 
interaction times. We further see that the stronger the interaction, the shorter the sub- 
Poissonian interval. Similar behaviour is available in a short interaction time if one turns 

(2) (2) 

his attention to g]^ (t), see Fig. 4b. However, in the large interaction times g^ (t) becomes 
stable in a sense of periodic oscillations. So, generally, we can see that for the Hamiltonian 
model (1.4), the modes lose their sub-Poissonian character and provide partially coherent, 
chaotic and superchaotic statistics for large interaction time. We should comment here that 
the behaviour of the 3rd mode in this model is different from that of the usual squeezed 
number states given in the literature earlier 2| even if they can provide a similar behaviour 
for the quadrature squeezing. This fact is clear by comparing the behaviour of the second- 



order correlation function of the 3rd mode with that of the usual number states l28ll where 
for the latter (with n ^ 0) g^'^^O) has smoothed behaviour and satisfies the inequality 
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(2) 

FIG. 4: Normalized normal second-order correlation function g- (t) against time for the modes 
which are initially in the number state |1,1,1) for a) the 1st mode for (Ai,A2,A3) = (0.6,0.7,0) 
(short-dashed curve), (0.6,0.7,0.7) (long-dashed curve), and (0.1,0.2,0.2) (solid curve); b) the 3rd 
mode for (Ai,A2,A3) = (0.6,0.7,0) (dashed curve), and (0.6,0.7,0.7) (solid curve). 

< ^^^^0) < 1-7. On the other hand, if one compares the short-dashed curves with the 
solid ones in Figs. 4a,b, it is easy to recognize that the two-mechanism interaction system 
ll| (A3 = 0) conserves the initial sub-Poissonian statistics of the interacting modes which 



can be periodically recovered with period t ^ 8 (for the 1st mode) and t ~ 16 (for the 3rd 
mode). Also similar approach is available for the model including competition between two 



parametric amplifiers and one parametric frequency converter 301]. 



The second quantity we are planning to study here is the violation of Cauchy-Schwarz 
inequality. This quantity shows the anticorrelation effects in the compound modes and can 
be observed in a two-photon interference experiment This quantity can be represented 



by the factor 41 1 



- 1. 



(3.10) 



{A]{t)Mt)Ai{t)A,{t)) 

The negative values for the quantity Ij,k{t) mean that the intermodal correlation is larger 
than the correlation between photons in the same mode |42|] and this indicates strong vio- 
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time 

FIG. 5: Evolution of the quantity Ij^kit) against the time for (nj,Xj) = (l,0.1),j = 1,2,3. The 
sohd and short-dashed curves represent the above quantity obtained in (1,2) and (1,3) modes, 
respectively. The long-dashed and star-centered curves are given for (Ai,A2,A3) = (0.1,0,0) and 
(0,0.1,0) , respectively, with fij = 1 . 

lation of the Cauchy-Schwarz inequality. 

The resulting expression for {A^j{t)Aj{t)A\{t)Ak{t)) is lengthy and will not be reproduced 
here. In Fig. 5 we have plotted the quantity Ij^kif) indicating the violation of Cauchy- 
Schwarz inequality between the jth mode and the fcth mode in dependence on time where 
{fij, \j) = (1, 0.1), j = 1, 2, 3. Further, the solid and short-dashed curves represent the quan- 
tity (3.10) obtained in (1,2) and (1,3) modes, respectively. For the sake of comparison the 
long-dashed curve is given for the single parametric amplifier, i.e. (Ai, A2, A3) = (0.1,0,0), 
for the same values of fij. Similarly the star-centered curve is given for the single frequency 
converter. In this figure the characters of both the Fock states and the structure of the 
Hamiltonian, (1.4), are reflected in the behaviour of the curves. To be more specific, the 
nonclassical negative values are dominant after switching on of the interaction. Also the 
behaviour of the h,2{t) is smooth and becomes stable for large values of interaction times, 
however, that of the h,3{t) undergoes amplification as well as oscillatory behaviour as a 
result of the energy exchange between these two modes in the system (cf. (1-4)). On the 
other hand, the comparison of the solid and long-dashed curves shows that the behaviour of 
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the /i,2(^) is close to that of the single parametric amplifier and for both /i,2(^) < for all 
times. This is in contrast with the behaviour of the two processes for the sum- squeezing, see 
Fig. 3. Nevertheless, this is not the case for the short-dashed (i.e. /i^3(t)) and star-centered 
curves where the former can exhibit anticorrelation only for a short time after switching on 
of the interaction, however, the latter evolves periodically with period t = 27r/A2 displaying 
always anticorrelation between modes. Further, similar behaviour can be seen regardless of 
the values of Xj. Finally, the behaviour of the l2,3{t) is similar to that of the /i,3(t). 



IV. QUASIPROBABILITY FUNCTIONS 

There are three types of quasiprobability functions: Wigner W-, Glauber P-, and Husimi 
Q-functions. These functions can be used as crucial to describe the nonclassical effects of 
the system, e.g. one can employ negative values of W^-function, stretching of Q-function and 
ligh singularities in P-function. Also, these functions are now accessible from measurements 



43| 



Indeed, the detailed statistics of the three coupled field modes can be obtained from 
photon-counting experiments. Most often we are interested in the quantum statistics of 
either one mode which determine the ensemble averages of the observables of this mode, 
or the composite statistics of the compound modes which refiect their mutually correlated 
properties, so that in the following we consider phase space distributions, in particular, 
VT-function for the single- and compound-modes when all modes are initially prepared in 
number states. To achieve this goal we calculate the s-parametrized joint characteristic 
function defined by 



C(3)(C,t,5) = Tr<^/3(0)exp 



3 



(4.1) 



where p(0) is the initial density matrix operator for the system under consideration, ( = 
(Ci)C2;C3)) s-iid s takes on values 1,0 and —1 corresponding to normally, symmetrically 
and antinormally ordered characteristic functions, respectively. The superscript 3 stands for 
three-mode case. When the modes are initially in number states \ni, n2, n^)', the symmetrical 
characteristic function can be calculated in a straightforward way as 

C^'\Cs = 0,t) = flexpi~\v,mLnA\vM")^ (4-2) 
where L„ is the Laguerre polynomial of order n and i]j{t) are given by 
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^i(^) = Ci/i W - C/2(t) - (293{t) + Qg,{t) + CMt) - Qh,{t), (4.3a) 
V2{t) = C29i{t) - C292{t) - CM) + CiUit) + Csh^it) - Qh,{t), (4.36) 
r/i(t) = Cshiit) - Qh2{t) + CiMt) - Cihit) - C295{t) + C296{t). (4.3c) 
We proceed by giving the definition of tlie s-parametrized joint quasiprobability functions 



as 



-I I l'c^'\Ct,s)llexpia,C;-a;Qd%, 
J J J • ^ 



(4.4) 



where a = (ai,a2,tt3)- When s = 1,0,-1, the formula (4.4) gives formally Glauber P- 
function, Wigner W^-function and Husimi Q-function, respectively. Now substituting from 



(4.2) into (4.4) and performing the integration we get the joint Wigner function as 



-1 



vni+n2+n3 



nexp(-2|6,(t)p)L„^,(4|6,(t)p 



(4.5) 



where €j{t) are 

ei(^) = "i/i(i) - «I/2(i) - a293it) + alg^it) + ash^^t) - alhe{t), (4.6a) 
e2(t) = alfi{t) - aif3{t) + a2fl'i(i) - a25'2(t) + a3^3(^) - alh^^t), (4.66) 
e3(t) = alfeit) - aif^{t) + algQ{t) - a29b{t) + a^hiit) - alh2{t). (4.6c) 
It is evident from (4.5) that the joint functions are strongly affected by the correlation 
between modes; this is reflected in the cross terms between different amplitudes of different 
modes. Such correlations have been used in a number of studies on nonclassical aspects of 
light including questions of violation of Bell inequalities 44 1. 

On the other hand, the ly-function of the single-mode can be obtained by means of 
integrating two times in the the joint ly-function (4.5) or by using the following relation 

iy«(«,t) = ^ j C^'\Q,t)expiaC;-a*Q)d%, (4.7) 
where C^^\Cj,t) is the single-mode symmetrical characteristic function which can be ob- 
tained from the joint symmetrical characteristic function (4.2). For instance, the character- 
istic function for mode Ai (say) can be obtained from (4.2) by simply setting C2 = Ca = 0- 
Now the single-mode Wigner function for the 1st mode when it is in the Fock states \ni) and 
the rest modes are in vacuum can be obtained from (4.7) together with (4.2) (after putting 
C2 

iy«(a,t) 



= in (4.2)), carrying out the integration, as 
2(-l)"i 



TT 



VA-m+it) 



exp 



2h+{t) 


ni 

E 

1=0 


\B-{t)] 


(ni-l) 




[A_{t)A^{t)\ 









(4.8) 



where 
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A^{t) = [hit) ± Mt)r + im ± u{t)r + m) ± ut)]', 

c±{t) = {a' + a*')[Mt)Mt) ± Mt)u{t) ± Mt)m], 



h±{t) = C^it) + D^it), 

2[A_(t) + B4t)][A4t)h4t) - B^{t)h+{t)] 



A^{t)B_{t)[A+{t)B.{t) - A^{t)B+{t)] 



^ 2[A+{t) + B^{t)][BAt)h+it) - A^{t)h_{t)] 

A4t)B4t)[A4t)B4t) - A4t)B^{t)] • 

The expressions for the 2nd and 3rd modes can be obtained from (4.8) using the transfor- 
mations (2.7). 

One can easily check when t — )■ that the expression (4.8) reduces to the well-known 



VT-function of Fock state \ni). In this case the following identity 45| 



ni 



5^L-(x)L-_,(y) = LZ-'^'^'i^ + y) (4.10) 
has to be used. 

Now we proceed by studying the evolution of the PF-function for the 1st mode when it is 
initially in the Fock state |1) and the other modes are in vacuum using the expression (4.8). 
In fact, the VT-function of the Fock state |1) is well known and it exhibits an inverted hole 
with more pronounced negative values related to the state showing maximum sub-Poissonian 
statistics. However, this behaviour can be completely washed out in our model at certain 
values of interaction parameters. This can be seen in Fig. 6a where one can observe the 
well-known shape for the l^-function of the thermal light, which is broader than that for 
vacuum and this demonstrates the Bose-Einstein statistics for thermal light exhibiting larger 
fluctuations than for coherent light. So the competition between different states in the model 
(power transfers) can cause that the output field exhibits classical nature of thermal light. 
In Fig. 6b, by increasing the values of the coupling constants, ly-function for squeezed 
thermal light dominates. This result is different from that given in section 3 since there the 



20 



(a) (b) 




(c) 



FIG. 6: VF-function against x and y for the 1st mode when it is in the Fock state |1) and the other 
modes are in vacuum for a) (t, Ai, A2, A3) = (|, 1.005, 1, 0.1) ; b) (t, Ai, A2, A3) = (|, 1.008, 1, 0.4) ; 
c) (t,Ai,A2,A3) = (tt, 0.5, 0.4, 0.4). 

initial thermal fields undergo squeezing when evolve in the present structure. It is important 
to mention that stretching here can be more pronounced for large interaction time. Finally, 
Fig. 6c displays an interesting fact that a statistical mixture for coherent superposition 
states, in principle, can also be realized where the basis of the two Gaussian bells dominates 
in the behaviour of 1^-function. This behaviour may be achieved also from the cat states 
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when they interact with a finite temperature heat bath or with the bath which consists 
of a gain medium in addition to the usual absorber j4^. For example, in the latter case this 
can be realized by choosing the gain appropriately so that one gets a purely diffusive motion 
of the field mode which leads to a double Gaussian structure with the missing oscillatory 
behaviour. Thus based on the fact that the nonlinear coupling between modes is present in 
the medium, the fields which are initially nonclassical cannot recover their properties during 
the evolution in the model for strong coupling and large interaction time. Nevertheless, 
they lose some nonclassical features, e.g. the I^-function loses its negative values, and is 
obtaining other ones, e.g. stretching of its form. 

0.20 
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FIG. 7: hoton-number distribution P{n^) for the 3rd mode when it is initially in the Fock state 
|1) and the other modes are in vacuum for t = 8 and (Ai, A2, A3) = (0.3, 0.3, 0.3). 



Finally, we would like to close this section by paying attention to the behaviour of the 
photon-number distribution when the modes are initially in Fock states. This will be done 
by means of the integral expression for the photon-number distribution in terms of Wigner 



function and Laguerre polynomials, which is given by 
^ ' ' W^(^)(a)exp(-2|anL„, 



Pin, 



(4.11) 



where j denotes the mode under consideration and W'^'^\a) is the single-mode Wigner 
function. We restrict ourselves to the photon-number distribution for the 3rd mode when 
it is initially in Fock state |1) and the other modes are in vacuum. After inserting the 
expression for VT-function of this case into (4.11) and performing the integral, the result has 
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been plotted in Fig. 7 for shown values. In this figure one can see the oscillatory behaviour 
typical for the photon-number distribution of squeezed states in the 3rd mode. This is clear 
in the pairwise nature of the oscillations, i.e. P{2n3) = 0, is positive integer. Indeed, for 
the squeezed states, i.e. squeezed vacuum, squeezed number and squeezed thermal states, 
these pairwise oscillations are explained as the result of the quadratic, or two-photon nature 
of the squeeze operator S{r) 



V. CONCLUSION 



In this article we have studied the statistical properties of three harmonic oscillators mu- 
tually interacting in the nonlinear crystal. The model is governed by interaction Hamiltonian 
including competition between two frequency converters and one parametric amplifier. Af- 
ter using the Heisenberg approach to the quantum statistics of interacting modes, we have 
investigated effects produced by the dynamics of the interaction as well as by the nonclassi- 
cal behaviour of the initial light modes. In our analysis we have considered that the modes 
are initially in the number states. We have proved that squeezed light can be generated 
in the standard sense or in the correlated quadratures (sum-squeezing). Furthermore, we 
have shown that this is still valid if the Fock states are replaced by thermalized fields. In 
this case the model operates as a microwave Josephson-junction parametric amplifier. We 
have also demonstrated that the initial sub-Poissonian statistics of the interacting modes are 
washed out, which leads to partially coherent, chaotic and superchaotic fields. Finally, we 
have shown that when the system is initially prepared in the Fock state |1,0,0), at certain 
values of interaction parameters, the 1st mode can evolve into thermal states or in squeezed 
thermal states. 
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